^ a SPECIFICATION 

^njKPg^ TITLE 

"METHOD FOR SURFACE-CONTOURING OF A THREE-DIMENSIONAL IMAGE" 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention concerns a method for surface contouring of a three- 
dimensional image of a subject. 
Description of the Prior Art 

For various applications (for example photographic analysis, industrial image 
processing, monitoring, traffic systems, media technology and animation), image 
processing and subject recognition of two-dimensional, and to a certain extent three- 
dimensional, images of subjects have progressed relatively far. The image data sets 
forming the images are, as a rule, relatively weakly disrupted and are acquired, for 
example, with optical sensors (see, for example, Inspect. MessTecAutomation, 
special publication 10/01, GIT Publishing Darmstadt, 2001, or Vision Systems 
Design, Vol. 6 #9, Pennwell, Sept. 2001). 

If the image data sets are relatively weakly disrupted, the surface contours of 
the images of the subjects can be relatively simply determined by means of known 
methods. With linear or nonlinear filtering, noise in the image data sets can be 
further reduced. In order to recognize textures of the images, pattern matching can 
be implemented. Borders of the images can be determined, for example, by 2D- or 
3D gradient formation. Techniques based on gradient formation are, for example, 
"active contours", specified in Chalana V., Linker D., "A Multiple Active Contour 
Model for Cardiac Boundary Detection on Echocardiographic Sequences", IEEE 
Trans. Med. Imaging 15, #3, June 1996, or what are known as "snakes", specified in 
Mclnerney T, Terzopoulos D: "T-snakes: Topology Adaptive Snakes", Medical Image 



Analysis 4, pages 73-91, 2000. Imaged subjects of the same brightness, color or 
texture can, on the other hand, be "filled out", for example by a technique known as 
region growing, and thus be completely described for a specific task (see Jendrysiak 
U., "Segmentierung von Schnittbildern". Spektrum der Wissenschaft, Dossier 1/1999: 
Perspektiven in der Medizintechnik 24-29, 1999). 

In medical technology, however, it may be the case that imaged healthy and 
diseased tissue exhibit a relatively small contrast difference, and image data sets are 
relatively severely noisy. This is true in particular for x-ray images (computed 
tomography, angiography, fluoroscopy), images produced with a magnetic 
resonance device or in nuclear medicine, and in particular images based on 
ultrasound (compare Marais P., Brady J., "Detecting the Brain Surface in Sparse 
MRI Using Boundary Models", Medical Image Analysis 4, 283-302, 2000, Chalana 
V., Linker D., "A Multiple Active Contour Model for Cardiac Boundary Detection on 
Echocardiographic Sequences", IEEE Trans. Med. Imaging 15, #3, June 1996, or 
Sakas G., "Dreidimensionale Bildrekonstruktion aus Ultraschall-Daten", Spektrum 
der Wissenschaft, Dossier 1/1999: Perspektiven in der Medizintechnik 18-24, 1999). 

In particular in a diagnostic finding for a patient, fine subject differences must 
be differentiated in the overall depth of the body or of the organs of the patient. 
Since the image data sets acquired with medical devices often are severely noisy, 
and acquired structures can be relatively complex, a determination of a contour in a 
two-dimensional image is already relatively computation-intensive. The 
determination of a contour of a three-dimensional image of a subject is even more 
difficult. Known methods therefore have refined, and often ad hoc adapted to, the 
present problem. An example is specified in Chalana V., Linker D., "A Multiple 
Active Contour Model for Cardiac Boundary Detection on Echocardiographic 
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Sequences", IEEE Trans. Med. Imaging 15, #3, June 1996. Another relatively 
computation-intensive method is specified in Jendrysiak U., "Segmentierung von 
Schnittbildem". Spektrum der Wissenschaft, Dossier 1/1999: Perspektiven in der 
Medizintechnik 24-29, 1999. 

SUMMARY OF THE INVENTION 
It is an object of the invention to provide a method with which a surface 
contouring of a three-dimensional image of a subject can be implemented with less 
calculation effort. 

This object is achieved in accordance with the invention by a method for 
surface contouring of a three-dimensional image of a subject including the steps of 
placing a first set of planes through the three-dimensional image, determining the 
contours that are imaged in each of the planes of the first set of planes and are 
associated with the surface contour of the three-dimensional image, and merging the 
contours determined in each plane of the first set of planes into a surface grid 
associated with the surface contour of the three-dimensional image. 

In accordance with the invention the surface of a three-dimensional image of 
a subject should be contoured or segmented. The three-dimensional image is stored 
in a volume data set, thus a data set of 3D blocks B(x,y,x) with KxLxM voxels 
(volume elements). The three-dimensional image, alternatively can be stored in the 
form of a data set that contains K slices of images with LxM pixels (picture elements) 
per image. 

For this, the planes of the first set of planes are placed through the three- 
dimensional image. The contour (that is associated with the surface contour of the 
three-dimensional image) is subsequently determined for each one of these planes. 
The determination of the contour associated with the individual planes corresponds 
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to the determination of contours of two-dimensional images. After this, the individual 
contours are merged into a surface grid associated with the surface contour of the 
three-dimensional image. 

In a preferred embodiment of the invention, the planes of the first set of 
planes intersect in a first straight line that, according to a particularly advantageous 
variant of the invention, proceeds substantially through the geometric center (center 
of gravity) of the three-dimensional image. This is particularly advantageous given a 
convex three-dimensional image. 

If, for example, the volume data set forming the three-dimensional image 
exhibits a relatively weak contrast, the volume data set can thus undergo a filtering, 
in particular for an automatic localization of the three-dimensional subject (according 
to a variant of the invention) for the determination of the geometric center of the 
three-dimensional image. The filtering in particular can be a median filtering, and the 
geometric center in particular can be automatically determined. The volume data set 
thus is significantly smoothed while preserving edges. Via such a significant 
smoothing, the detailing of the boundaries is lost, however the geometric center of 
the three-dimensional image can be relatively precisely determined. 

If the subject to imaged is a part of a living organism, and the three- 
dimensional image is acquired by means of a medical technology device, the 
geometric center of the three-dimensional image can be determined during a 
navigation-guided treatment of the living organism, according to a further 
embodiment of the invention. The geometric center can be determined, for example, 
during the navigation-guided treatment by the position of the focus of a navigated 
medical instrument. 
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A medical device suitable for the navigation is, among other things, specified 
in German OS 199 51 502. For example, a volume data set of the region of interest 
of the body of a patient is generated with an x-ray device. In particular during a 
minimally invasive procedure, for example, a surgeon inserts the medical instrument 
into the body of the patient. A position sensor of the navigation system is arranged 
at the medical instrument, such that the navigation system determines the position 
(meaning the location and orientation of the medical instrument) relative to the body 
of the patient. Based on the position determination, thus based on a determination 
of the spatial coordinates of the medical instrument, an image of the medical 
instrument can be blended into the image associated with the volume data set (that 
is, for example, shown on a monitor). 

An automatic recognition of the approximate position of the medical 
instrument also is useful in nuclear medicine. A volume data set acquired by PET 
exhibits a relatively good contrast, however it has a relatively poor resolution. An 
alternative is to manually designate the geometric center with a mouse click in 
another representation, for example in an x-ray image produced using contrast 
means in two projections, or in a CT-MPR multi-planar reformatting representation. 
This is nothing more than the reassembly of the volume data set in another 
orientation that, for example, originally horizontal slices are re-oriented (orthogonal 
MPRs: horizontal, sagittal, coronal; free MPR: arbitrary slanted secondary slices). 
This works isotropically with new multi-line spiral CTs or MR devices, for example 
from 512 CT slices with 51 2 2 pixels per slice). The background for this is that 
various views of the image of the subject and its surroundings are shown coupled in 
a defined position and orientation. If the three-dimensional view is shifted and/or 
rotated, the MPR representation oriented thereto moves with it. The three- 
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dimensional shifting and rotation are continued until the image of interest is centrally 
intersected by the MPR slice. Then the middle of this planar representation of the 
image is clicked; the subject center is thus defined in three dimensions. 

By means of the inventive method, the dimensionality of the volume 
processing (3D image processing) is reduced by at least one stage. The curve of a 
3D subject density D of the three-dimensional image along a coordinate axis x of the 
space can be graphically represented as a two-dimensional function D1(x). The 
three-dimensional image in space therefore can be correspondingly described with 
a function D3(x,y,z) whose values (characteristic variables) are plotted in the fourth 
dimension. Thus the space spanned by x, y and z can be projected such that only 
the problem similar to the "one-dimensional wandering" (see Tittmann P., "Irrfahrten 
in Graphen", Memo Hochschule Mittweida, Sept. 1998) or, respectively, an 
optimization along an abscissa x' alone (instead of (x, y, z)) is to be solved. 

In order to achieve this, the planes of the first set of planes are placed through 
the three-dimensional image. In the general case, the planes of the first set of 
planes should if at all possible not redundantly, but completely, cover all voxels of 
the three-dimensional image in various orientations. The planes of the first set of 
planes preferably intersect in the first straight line that, at least substantially, 
proceeds through the geometric center of the three-dimensional image. With regard 
to a fixed normal axis, the individual planes of the first set of planes can exhibit, for 
example, equidistant angles ch through a m , when the first set of planes comprises m 
planes. By the transition in a rotation geometry, redundancies arise in the direction 
of the rotation point, and lesser scanning results in the opposite direction. However, 
scanning primarily needs to be observed in the region of the expected surface 
contour. By the "new acquisition" of the three-dimensional image in planes rotating 
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opposite to one another that intersect the three-dimensional image, a first stage of a 
dimension degradation is achieved, since the data sets associated with the planes 
are two-dimensional data sets and, for example, are described with Cartesian 
coordinates. If the volume data set is described with Cartesian coordinates (x,y,z), 
each data set associated with the planes of the first set of planes thus can be 
described with two-dimensional Cartesian coordinates (Xj",yj"), whereby j=1 ... m, 
when the first set of planes contains m planes. Subsequently, the data sets 
associated with the individual planes of the first set of planes are analyzed in order to 
determine contours that are imaged in each of the planes of the first set of planes 
and are associated with the surface contour of the three-dimensional image. 

In each of these planes (Xj",yj"), a closed contour is expected. The contour 
K," of the j-th plane of the first set of planes has the points K/' = (xj/'.yjj"), with i = 1 
... n. In order to determined each Kjj", a two-dimensional search is required in 
( x j,i">yj/') that is known in general to those skilled in the art. 

In order to reduce the calculation effort, however, in an embodiment of the 
invention, when the image contents of each plane of the first set of planes are 
specified with Cartesian coordinates, the following steps are implemented for each 
plane for the determination of the contours of each individual planes. A coordinate 
transformation is implemented according to polar coordinates with regard to the 
geometric center of the three-dimensional image, and therewith unwinding of the 
contour, and the contour is determined in the transformed plane. 

The contour in the transformed plane is thereby determined according to a 
preferred variant of the invention, in which a contour continuation ensues in the 
direction of the angular coordinate (<t>) with an optimization according to minimal 
changes. 
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Thus a transformation according to polar coordinates is applied with regard to 
the geometric center and the first straight line. A "contour field" K'(x',ct>) with fixed <3> 
thus results that covers 360°. The radial x' is designated as r (r=x') in the following, 
corresponding to the polar coordinates. With relatively little effort, only the one- 
dimensional problem remains, to determine the radius r at each angle O as a single 
quantity. The fact that this is actually a rather "one-dimensional" problem can be 
illustrated by plotting this function on a Cartesian system, for example with cp in the 
abscissa and r in the ordinate. Some boundary conditions now result for r that can 
be utilized in order to find a qualitatively good contour: 

The contour must have the same r-coordinate for Oi and O n (closed 

contour). 

The contour is continuous and has a predominantly vertical gradient. 
The r-offset from line to line is relatively small. 

A typical "edge profile" in the brightness profile can be found by line 
comparison. 

Possible compensation of signal noises in the line comparison. 
For the transformation according to polar coordinates, the angular increment 
should preferably be selected such that adjacent points (Oj, r) and r) in the 

region of the contour have approximately the same separation of a voxel edge 
length. In place of the transformation according to polar coordinates, a 
transformation with regard to a "model contour" can also ensue when a similarity to 
such a model is to be expected. The abscissa is then characterized by the 
continuous indices i1 through in of the contour points of the model; the perpendicular 
separations of the model contour are plotted on the ordinate. Here the 
transformation according to polar coordinates is a rescanning. 
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In the matrix resulting from the rescanning, various features for an object 
edge can now be determined or accentuated by suitable signal processing, as well 
as being brought into context. For example, a search method can be applied that 
analyzes and combines all reasonable path branchings. An example of a realization 
is the dynamic optimization according to R. Bellman, "Dynamic programming and 
stochastic control processes", Information and Control, 1(3), pages 228-239, 
September 1958. The result is a continuous contour, for example in the polar 
coordinate representation of O-i through 4> n (or, respectively, for h through i n ). This 
contour can be transformed back into the original Cartesian coordinates. The 
dynamic optimization can combine a number of feature values (filter responses) into 
a suitable cost function, and, according to the forward accumulation of costs, 
includes a backtracking as a substantial element. A plane of the first set of planes is 
thereby contoured. 

The remaining planes of the first set of planes are correspondingly analyzed. 
Already-determined contours of previously analyzed planes can be used in order to 
also ensure the consistency of the three-dimensional surface contour of the three- 
dimensional image in the rotation direction described by a (usable prior knowledge). 
Likewise, all planes of the first set of planes can be analyzed independently. With 
one or more successive processing stages with which the individual contours of the 
planes of the first set of planes are compared, the obtained surface grid can be 
refined. 

At the end of the processing of all planes of the first set of planes, the surface 
contour of the three-dimensional image initially exists in coordinates (r, a, O). In this 
"cuboid", a processing of the 3D connection [correlation; interrelation] is possible, for 
example as a filtering of the r-values with a 2D operator. 
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The transformation back into the original Cartesian coordinates yields the 
surface grid of the three-dimensional image. The surface grid is comparable to a 
net-like cage. The surface grid now can be improved and refined in the original 
coordinates. The voxels within the surface grid can subsequently be marked as 
belonging to the three-dimensional image. 

The inventive method in the context of a spherical geometry has the 
advantage that it possesses no preferred direction and that, considered statistically, 
in practice occurring surface gradients are handled better on average than in a 
method that is based on a Cartesian grid, since such a surface lies more parallel to a 
closed 3D contour. Because the inventive method is, geometrically, relatively well- 
adapted to the problem, the number of the planes of the first set of planes can be 
smaller by a multiple than when, in the region of the surface, the surface would be 
sampled in the voxel separation; this means the "net lines" can have a separation of, 
for example, ten or more voxel edges from one another while still describing the 
overall contour with sufficient precision. 

In particular for an improved determination of the surface contouring of the 
three-dimensional image, or when unfavorable output data are present, according to 
an embodiment of the invention the additional steps are implemented. 

A second set of planes is placed through the three-dimensional image, with 
the planes of the second set of planes intersecting in a second straight line. The 
contours that are imaged are determined in each of the planes of the second set of 
planes and are associated with the surface contour of the three-dimensional image. 
Together with the contours determined in each plane of the first set of planes, the 
contours determined in each plane of the second set of planes are merged into the 
surface grid associated with the surface contour of the three-dimensional image. 
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By re-extraction oriented on net lines (scan lines proceeding at right angles to 
one another, with orientation from the inside out) determined by means of the planes 
of the first set of planes and associated with the surface grid, the contouring also 
works with less convex, irregular three-dimensional images. 

A benefit of reliability and precision can be achieved when, according to an 
embodiment of the invention, the first straight line is aligned at a right angle to the 
second straight line, and in particular the first straight line and the second straight 
line intersect in the geometric center of the three-dimensional image. 

The precise surface contouring can be used, for example, to determine the 
three-dimensional image precisely for a functional evaluation. For therapeutic or 
surgical procedures, the precise shape in relation to the vicinity (for example 
vessels), and also the shape change in the course of time, is very important. Further 
applications are possible in plastic surgery or in prosthetics, for the recognition of the 
subject to be replaced or of cavities for an implant. 

A simplified calculation of the necessary calculation operations illustrates 
some of the basic advantages of the inventive method. It is assumed that the 
volume data set exists in the form of 50 computed tomography slices of the matrix 
512x512. The three-dimensional image is formed by a sub-cuboid of 40x200x200 
voxels. 

Given a surface contouring of the three-dimensional image according to the 
prior art, starting from a start point its continuation on the surface is sought. The 
surface of the entire sub-cuboid is 112000 pixels. For the gradient formation, a 
factor of 1 1 is applied with regard to the surface (coupler field 3 3 =*"26"), and applied 
with a factor of "1" for the 27 covered voxels. This effectively corresponds to five 
surface search-stripes inside and outside of the surface. Starting from each point, in 
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principle 17 directions in a preferred direction and laterally are examined (start) that 
then reduce inside a surface. In a more sensible manner, the search is empirically 
delimited as specified, for example, in Buckley K., Vaddiraju A., Perry R., "A New 
Pruning/Merging Algorithm for MHT Multitarget Tracking", Radar-200, May 2000. 
Furthermore, it is assumed that it is sufficient to on average examine six alternatives 
with an average depth of three in order to differentiate a reasonable continuation 
from a meander or, respectively, to land on an already-classified point (multiplier 6 3 = 
216). 

Given a surface contouring according to the already specified method, and 
not under consideration of "backtrackings", approximately 11x112000x26 = 
33264000 gradient formations are required, and approximately 216x112000 = 
24192000 calculation operations for the search for maximum gradient sum. Thus 
approximately 57000000 calculation operations result. 

For the approximate calculation of necessary calculation operations for the 
surface contouring of the same three-dimensional image according to the inventive 
method, it is assumed that the first set of planes comprises 20 planes. After the 
transformation into polar coordinates, the image information of the transformed 
planes are respectively sampled in the 2° raster, whereby the corresponding CT 
images are reconstructed in the typical raster 512x512: 180x256 = 46000. The 
gradient formation and contour search are first implemented in the planes. 
Moreover, a noise-suppressing folding (fold-over) operation is implemented with a 
coupler field of seven. A typical density profile in the region of the subject edge is 
generated with seven supporting locations and shifted over the entire radius. In 
each shifting position, the covariance is calculated as an iteration value. Given the 
search of the contours, respectively a "jump" by one pixel from line to line is 
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provided; there are thus three alternatives. For post-processing of the connection 
between the planes, a 3x3 coupler field is used. It thus results in approximately 
20x46000 = 920000 operations for the gradient formation, and approximately 
3x46000 = 138000 calculation operations for the search for a maximum gradient 
sum. Approximately 20x180x9 = 32400 calculation operations result for the 
backtracking, which here corresponds to a testing of the plane connection. Thus 
approximately 1090400 calculation operations result. 

Based on the inventive method with the cited assumptions, a reduction of the 
effort in calculation operations results by approximately the factor of 50. 

DESCRIPTION OF THE DRAWINGS 

Figure 1 shows a computed tomography apparatus operable in accordance 
with the invention. 

Figure 2 shows a set of planes whose planes intersect with a three- 
dimensional image acquired with the computed tomography apparatus shown in 
Figure 1. 

Figure 3 shows, in Cartesian coordinates, image information of one of the 
planes shown in Figure 2. 

Figure 4 shows, in polar coordinates, image information of the plane shown in 
Figure 3. 

Figure 5 is an illustration of a determination of a contour shown in Figure 4. 
DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Figure 1 schematically shows a computed tomography apparatus with an x- 
ray source, from which is emitted a pyramid-shaped x-ray radiation beam 2 (the 
edge rays of which are shown as dash-dot lines in Figure 1) that proceeds through 
an examination subject (for example a patient 3) and strikes a radiation detector 4. 
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In the exemplary embodiment, the x-ray source 1 and the radiation detector 4 are 
arranged opposite one another on an annular gantry 5. With regard to a system axis 
6 that runs through the middle point of the annular gantry 5, the gantry 5 can be 
rotated on a mounting device (not shown in Figure 1) (arrow a). 

In the exemplary embodiment, the patient 3 lies on a table 7, transparent to x- 
ray radiation, which is movably positioned along the system axis 6 by means of a 
carrier device (likewise not shown in Figure 1 ) (arrow b). 

The x-ray source 1 and the radiation detector 4 thus form a measurement 
system that, with regard to the system axis 6, is rotatable and movable along the 
system axis 6 relative to the patient 3, such that the patient 3 can be irradiated at 
various projection angles and various positions with regard to the system axis 6. 
From the ensuing output signals of the radiation detector 4, a data acquisition 
system 9 thus forms measurement values that are supplied to a computer 11 that 
calculates, by means of method known to those skilled in the art, an image of the 
patient 3 that can be reproduced on a monitor 12 connected with the computer 11. 
In the exemplary embodiment, the data acquisition system 9 is connected with an 
electrical line 8 that has , for example, (in a manner not shown) a slip or wiper ring 
system, or a wireless transmission link, with the radiation detector 4, and with an 
electrical line 14 with the computer 1 1 . 

The computed tomography apparatus device shown in Figure 1 can be used 
both for sequential scanning and for spiral scanning. 

In sequential scanning, a slice-by-slice scanning of the patient 3 ensues. The 
x-ray radiation source 1 and the radiation detector 4 are thereby rotated around the 
patient 3 with regard to the system axis 6, and the measurement system comprising 
the x-ray source 1 and the radiation detector 4 acquires a plurality of projections in 
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order to scan a two-dimensional slice of the patient 3. A sectional image 
representing the scanned slice is reconstructed from the thereby-acquired 
measurement values. The patient 3 is moved along the system axis 6 between the 
scanning of successive slices. This event repeats until all slices of interest are 
acquired. 

During spiral scanning, the measurement system formed by the x-ray source 
1 and the radiation detector 4 rotates with regard to the system axis 6, and the table 
7 continually moves in the direction of the arrow b, meaning the measurement 
system formed by the x-ray system 1 and the radiation detector 4 continually moves 
relative to the patient 3 on a spiral track c until the region of interest of the patient 3 
is completely acquired. A volume data set is thereby generated that, in the case of 
the present exemplary embodiment, is coded according to the DICOM standard 
typical in medical technology. 

In the exemplary embodiment, the surface of a three-dimensional image of a 
subject (acquired with the volume data set) to be contoured is part of the region of 
interest of the patient 3. In this case, the subject is a tissue variation of the patient 3. 
The three-dimensional image 21 is schematically shown in Figure 2 and can be 
viewed with the monitor 12. 

In order to contour the three-dimensional image 21, a first set of planes that 
(in the exemplary embodiment) contains 20 planes Ei through E 2 o is placed through 
the three-dimensional image 21. For clarity, in Figure 2 only the planes E1 through 
E5 are shown. The planes E1 through E20 intersect in a first straight line G that, in 
the exemplary embodiment, proceeds through the geometric center S of the three- 
dimensional image 21. The geometric center is either manually, approximately 
marked in the three-dimensional image 21 by an operator (not shown in Figure 1) by 
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means of a mark (not shown in Figure 2) that can be mixed into the field of view of 
the monitor 12 and that can be moved with a computer mouse 13 connected with the 
computer 11, or the computer 11 is (in the exemplary embodiment) configured such 
that it automatically determines the approximate geometric center S of the three- 
dimensional image 21, and subsequently automatically determines the planes Ei 
through E 2 o- Moreover, the operator in the exemplary embodiment can subject the 
volume data set to, for example, a median filtering with the computer 11. The 
volume data set is thereby significantly smoothed while preserving edges, in order to 
improve the determination of the geometric center. 

In the exemplary embodiment, the planes E1 through E20 rotate opposite one 
another, respectively by an angle Aa = 9°, and the first straight line proceeds parallel 
to the x-axis of the DICOM coordinate system already defined for the volume data 
set. If, for example, the three-dimensional image 21 has a quarter of the overall size 
of the isotropic volume data set in the exemplary embodiment, the maximum point 
separation at the surface of the three-dimensional image 21 thus corresponds to 
approximately 20 pixels. 

In each of the planes E1 through E20, slices of the three-dimensional image 
21 are imaged. Such a section 31, which in the exemplary embodiment is fixed 
through the plane E1, is schematically shown in Figure 3. Moreover, the image 
contents of the individual planes Ei through E 2 o are described in Cartesian 
coordinates. The image contents of the plane E 1 are described with the Cartesian 
coordinates x" and y" in the exemplary embodiment. 

In order to obtain the surface contour of the three-dimensional image 21, the 
slices shown in the planes E-! through E 2 o, thus also the slice 31 of the plane Ei, are 
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subsequently analyzed, in that the contours of the individual slices are determined. 
The contour of the slice 31 is provided with the reference character 32 in Figure 3. 

For an efficient determination of the contours, in the exemplary embodiment 
the image information of each plane through E 2 o is radially newly arranged, by 
being transformed with regard to the geometric center S according to polar 
coordinates (r, O). The result of such a coordinate transformation is shown as an 
example in Figure 4, for the transformation of the image contents in particular of the 
slice 31 of the plane E1 shown in Figure 3. The transformed contour of the slice 31 
is provided with the reference character 32' in Figure 4. 

The transformed contour 32' is subsequently determined. A limitation path 
(definition path) is thus sought that corresponds to the contour 32' transformed into 
polar coordinates. For this, a direction-focuses filtering first ensues that can even be 
one-dimensional in the r-direction, and that can operate so as to both be sensitive to 
edges and to suppress noise. A one-dimensional search for the limitation path, thus 
for the transformed contour 32', now is undertaken using filtering according to the 
principle of dynamic optimization, as, for example, specified in R. Bellman, "Dynamic 
programming and stochastic control processes", Information and Control, 1(3), 
pages 228-239, September 1958. Thus the r-coordinates corresponding to the 
angular coordinates <t> between 0° and 360° are sought, with the boundary condition 
that the r-coordinate for O = 0° should be equal to the r-coordinate for O = 360°. In 
the exemplary embodiment, these r-coordinates are determined by, after a relatively 
significant smoothing, superimposing the first and last five lines, and subsequently 
the edge positions are determined in the smoothed line data. This search is shown 
as an example in Figure 5, wherein the transformed contour is again provided with 
the reference character 32\ The lines 40 in Figure 5 are shown thinner than the 
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transformed contour 32' and represent path alternatives found during the 
determination that, however, were again discarded. 

The contours (transformed into polar coordinates) of the remaining planes E 2 
through E 20 are similarly determined and subsequently considered in context. This 
involves examining whether the continuity of the brightness or of the gradients 
relevant in this direction is plausible. 

The found contours (transformed into polar coordinates) of each plane Ei 
through E 2 o are subsequently transformed back into Cartesian coordinates, and after 
this are merged into a surface grid. The surface grid thereby corresponds to the 
surface contour of the three-dimensional image 21 . 

In the exemplary embodiment, the surface grid is now improved and refined in 
the original three-dimensional Cartesian coordinates. Using the surface grid, the 
surface contour of the three-dimensional image 21 to be determined is interpolated 
up to the point separation. The voxels within the surface grid are subsequently 
marked as pertaining to the three-dimensional image 21 . 

The volume data set is shown as a three-dimensional data block B'(r, O, a) by 
the placement of the planes Ei through E 2 o through the three-dimensional image 21 
and via the subsequent coordinate transformation according to polar coordinates. 
The determination of the transformed contours (such as, for example, the 
transformed contour 31') is less calculation-intensive for approximately convex 
subjects. A contiguous "mountain surface" is sought via O and a. The "mountain 
surface" locally exhibits relatively minor oscillations in the r-direction. 

The number of the voxels analyzed for the determination of the contours is 
also reduced with regard to a three-dimensional surface contouring according to the 
prior art. If the volume data set contains, for example, 50 CT slices with an 
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approximate dimension of 128 pixels of the three-dimensional image whose surface 
contour should be determined (thus a processing field approximately 50x256x256 
large), approximately a million voxels (20 planes, polar coordinate raster 2°, r max = 
256) are analyzed. According to the prior art, in contrast to this approximately 3 
million voxels are analyzed. 

In order to obtain the improved surface grid in the exemplary embodiment, in 
addition to the planes Ei through E 2 o of the first set of planes, alternatively at least 
one further, second set of planes (not shown in the Figures for reasons of clarity) can 
be placed. The planes of the second set of planes intersect in a second straight line 
that likewise proceeds through the geometric center S and is aligned at a right angle 
to the first straight line G. The slices shown with these planes are again transformed 
into polar coordinates with regard to the geometric center S. The transformed 
contours that are associated with the surface contour of the three-dimensional image 
21 are subsequently determined again, analogously as above. The determined 
transformed contours are subsequently again transformed back, and thereafter are 
merged into the surface grid. 

It is also an advantage that, via the coordinate transformation into polar 
coordinates, the determination of the transformed contours in the rotation direction a 
(azimuthal) and in directions within the planes (ct>, orbital) can substantially be three- 
dimensionally implemented as a two-dimensional filtering, namely as a smoothing of 
r(<t>,a), instead of in the original coordinates. 

The volume data set with which the three-dimensional image 21 is produced 
is, in the exemplary embodiment, produced with a computed tomography apparatus. 
The volume data set alternatively can be produced with other medical technology 
devices, such as in particular with an x-ray device, a magnetic resonance device, by 
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PET, an ultrasound device, or by means of nuclear medicine. Imaging devices or 
technical devices outside of medical technology also can be used. 

The three-dimensional image 21 is an image of a part of the patient 3, 
however, the surface of a non-human three-dimensional image also can be 
contoured by means of the inventive method. 

The first straight line G need not necessarily proceed through the geometric 
center S of the three-dimensional image 21. The planes Ei through E 2 o also need 
not necessarily intersect in the straight line G. 

Although modifications and changes may be suggested by those skilled in the 
art, it is the intention of the inventor to embody within the patent warranted hereon all 
changes and modifications as reasonably and properly come within the scope of his 
contribution to the art. 
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